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ABSTRACT 

We study the generalized Kompaneets equation (kinetic equation for the photon distribution 
function taking into account the Compton scattering by electrons) using a relativistically 
covariant formalism. Using the generalized Kompaneets equation, we derive an analytic 
expression for the Sunyaev-Zel'dovich effect, which takes into account up to 0(9^) terms, where 
Oe = kBTe/mc^ is the relativistic expansion parameter and is the electron temperature. We 
carefully study the applicable region of the obtained analytic expression by comparing with the 
result of the direct numerical integration. We conclude that the present analytic expression can 
be reliably applied to the calculation of the Sunyaev-Zel'dovich effect for fes^e < 15keV, which 
is the highest electron temperature in the presently known clusters of galaxies. Therefore, the 
present analytic expression can be applied to all known clusters of galaxies. 

Subject headings: cosmology: theory Hubble constant cosmic microwave background 
radiation — galaxies: clusters: Sunyaev-Zel'dovich effect — plasmas: Compton scattering 
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1. INTRODUCTION 

Compton scattering of the cosmic microwave background (CMB) radiation by hot intracluster gas — 
the Sunyaev-Zel'dovich effect — (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972, 1980) provides 
a useful method to measure the Hubble constant Hq (Gunn 1978; Silk & White 1978; Birkinshaw 1979; 
Cavaliere, Danese, & De Zotti 1979; Birkinshaw, Hughes, & Arnaud 1991; Birkinshaw & Hughes 1994; 
Myers et al. 1995; Herbig et al. 1995; Jones 1995; Markevitch et al. 1996; Holzapfel et al. 1997). The 
Sunyaev-Zel'dovich formula has been derived from a kinetic equation for the photon distribution function 
taking into account the Compton scattering by electrons — the Kompaneets equation — (Kompaneets 
1957; Weymann 1965). The Kompaneets equation has been derived with the non-relativistic approximation 
for the electron. However, the electrons in the clusters of galaxies are extremely hot, fcsTg = 5 15keV 
(Arnaud et al. 1994; Markevitch et al. 1994; Markevitch et al. 1996; Holzapfel et al. 1997). 

Recently attempts have been made to include the relativistic corrections in the Sunyaev-Zel'dovich 
effect (Rephaeli 1995; Rephaeli & Yankovitch 1997). However, it appears that the calculations have not 
been carried out in a manifestly covariant form. For example, equation (4) in Rephaeli (1995), which 
comes from Chandrasekhar (1950), is a non-relativistic formula. Since the extension of the Kompaneets 
equation to the relativistic regime is extremely important in view of many recent measurements of the 
Hubble constant Hq with the use of the Sunyaev-Zel'dovich effect, we will solve the kinetic equation for the 
photon distribution function in a manifestly covariant form taking into account the Compton scattering by 
electrons. 

Very recently a generalized Kompaneets equation has been derived by two groups (Stcbbins (1997); 
Challinor & Lasenby (1997)). By using the generalized Kompaneets equation, analytic expressions for the 
Sunyaev-Zel'dovich effect have been derived as a power series of 9^ = kBTe/mc? , where T,. and m are the 
electron temperature and the electron mass, respectively. It has been shown that the results obtained by a 
power series expansion agree with the previous numerical calculations by Rephaeli (1995) and Rephaeli & 
Yankovitch (1997). Analytic expressions are compact and extremely useful to study the Sunyaev-Zel'dovich 
effect. On the other hand, it has been pointed out by Challinor & Lasenby (1997) that the convergence 
of the power series expansion in is slow. They suspect that the expansion is an asymptotic expansion. 
Therefore, it is extremely important to study the extent of the 6e region where the analytic expressions 
can be applied. This is the main subject in the present paper. Following the approach of Challinor & 
Lasenby (1997), we will derive analytic expressions for the intensity change of the photon spectrum. In the 
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derivation we will take into account the relativistic terms up to the fifth order in Og. In order to examine 
the accuracy of the analytic expressions derived from the power series expansion, we will directly integrate 
the Boltzmann equation numerically. Comparing these results we will carefully study the valid region for 
the electron temperature and for the photon angular frequency oj where the present analytic expressions 
can be reliably applied. 

The present paper is organized as follows. Extension of the Kompaneets equation to the relativistic 
regime will be treated and an analytic expression for the intensity change of the photon spectrum will 
be derived in § 2. In § 3 we will study the accuracy of the analj^tic expressions by comparing with the 
numerical results. Concluding remarks will be given in § 4. 



2. GENERALIZED KOMPANEETS EQUATION 

In this section we will extend the Kompaneets equation to the relativistic regime. We will formulate 
the kinetic equation for the photon distribution function using a relativistically covariant formalism 
(Berestetskii, Lifshitz, & Pitaevskii 1982; Buchler & Yueh (1976). As a reference system, we choose the 
system which is fixed to the center of mass of the cluster of galaxies. This choice of the reference system 
affords us to carry out all the calculations in the most straightforward way. We will use the invariant 
amplitude for the Compton scattering as given by Berestetskii, Lifshitz, & Pitaevskii (1982) and by Buchler 
& Yueh (1976). 

The time evolution of the photon distribution function n{u}) is written as 

^ = -2 / -0^d'p'd'k' W {n{uj)[l + n{u:')]f{E) - n{J)[l + n(a;)]/(£')} , (2.1) 

2oj(jj'EE' ' ^ ' ' 

K = -2{p ■ k) = -2ujE ^1 - i^cosa^ , (2.4) 

k' = 2{p ■ k') = 2u'E (l - -^cosa'^ . (2.5) 

In the above W is the transition probability corresponding to the Compton scattering. The four-momenta 
of the initial electron and photon are p = {E,p) and k — (cj, fc), respectively. The four- momenta of the final 
electron and photon are p' = {E',p') and k' = (u)', k'), respectively. The angles a and a' are the angles 
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between p and k, and between p and k' , respectively. Throughout this paper, we use the natural unit 
h = c = 1 unit, unless otherwise stated explicitly. 

By ignoring the degeneracy effects, we have the relativistic Maxwellian distribution for electrons with 
temperature Te as follows 

-1 



m = 



^{{E-m)-{^-m)}/kBTe ^ ^ 
-{K-{n-m)}/kB'l\. 



(2.6) 



where K = {E — m) is the kinetic energy of the initial electron, and (/U — m) is the non-relativistic chemical 
potential of the electron. We now introduce the quantities 



X = 



Ax 



ksTe 



ksTe 



Substituting eqs. (2.6) - (2.8) into eq. (2.1), we obtain 



dn{u)) 



= -2 



d^p'ifk'W f{E) {[1 + n{Lo')]n{u) - [1 + n{io)]n{Lo')e^'=) . 



(2.7) 
(2.8) 

(2.9) 



dt ~ J (27r)3 

Eq. (2.9) is our basic equation. One can numerically integrate this equation directly. We will perform this 
integration in §3. 

Following Challinor & Lasenby (1997), we expand eq. (2.9) in powers of Aa; by assuming Aa; <C 1. We 
obtain the Fokker-Planck expansion 



dn{ui) 
dt 



= 2 

+ 2 

+ 2 

+ 2 

+ • • 



On 



+ n{l + n) 



h 

dn 



dx 

^+2(l + n)|^+n(l + n) 

d^n „/, sd'^n .dn 
^+3(l+n)— +3(l + n)-+n(l + n) 

d^n .d^n .d^n .dn 

^+4(l+n)^+6(l+n)^+4(l+n)-+n(l + n) 



(2.10) 



where 



k\l 



_fp_ 



d-'p'd-'k'W f{E) {Axf . 



(2.11) 



Analytic integration of eq. (2.11) is not possible except for doing power series expansions of the integrand. 
Technically speaking, there are several choices for the expansion parameter of the integrand of eq. (2.11). 
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They are, for example, p/m, K/m = {E/m — 1), and v = p/E. It is important to note that obtained 
analytic expressions of Ik after the integration do not depend on the choice of the expansion parameter. It is 
also extremely important to note that the expansions in terms of these variables are asymptotic expansions 
in 7fe. Therefore, not only the convergence is very slow but also the accuracy of the analytic expressions 
have to be carefully examined for the parameter region considered. This is one of our main subjects in the 
present paper. 

Challinor & Lasenby (1997) carried out a calculation up to 0{9l) terms. We will carry out a calculation 
up to O{0l) terms in the present paper. The calculation of Ik has been performed with a symbolic 
manipulation computer algebra package Mathematica. We obtain 

1 f / 47 21 

h = -(TTNeOeX i4.-x + eAlO-—x + —X' 

,2 /15 1023 567 2 147 , 

15 2505 9891 ^ 9551 , 1616 , 
2 8 10 20 35 

^,,'135 30375 177849 o 472349 , 63456 , 940 , , , , , 

+ ^e(^-^X+^.^-^.3^__,4__^5,^, (2.12) 



1 n 2 \ . ^ /^47 63 7 2 

= -^-TN.(^ex'^l + e^[^---x+-x^ 

^2 /1023 1302 161 2 22 



+ ol 



+ o\ 



\ o o ^ "J / 

/2505 10647 38057 , 2829 , 682 A 

j 

/30375 _ 187173 1701803 _ 44769 61512 _ 510 A 1 

V 128 20 80 4 35 7 y J ' ^ ' 



1 , r /42 7 

-aTNeOeX^ < 9e{—--X 



2 I V 5 5 

,2 /'868 658 88 2 H ^ 
^•^ ' -X H x^ x'^ 







5 5 30 



/ 7098 14253 8084 „ 3503 , 64 , 

+ '"[-5 +21^ 

r 62391 614727^ ^ _ 123083^3 ^ 14^,4 _ 3^,. , , ^ (^.14) 



I- 



10 20 16 21 21 



74 = larNJex'^ , 

2 I 10 
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/« = 



^2 /329 11 2 

^, / 14253 9297 7781 o 320 , 16 . 
+ ,4 ^614727 _ 12^89^ ^ 239393^2 _ 1^,3 ^ 12676^, _ H ^ ^ ^ ^^^^^ 



1 /44 11 



2 I \ 5 10 



,3 / 18594 36177 192 ^ 64 3 



+ 0t X H a;^ x' 

" ' 35 140 7 105 

^4/ 124389 1067109 2 5032 3 66 . H r; , , 



^3 / 12059 64 32 ^\ 



/ 355703 8284 6688 2 3 H 4 1 1 

a; + —— a;^ - 22a;3 + — n- , (2.17) 



V 80 3 15 42 

^UTN^Oex' I el 

( 



/7 = ;:(TTN,eex' ^eu^-^x 



,,^16568 30064^^ 176^2. n 3 ^^^^^^ 



21 105 7 21 



1 ^ 8 f 16 .3 .4 / 7516 99 11 2M 



I, = 1 TVA^^I^^f^-lixlk (2.20) 



'10 



2 V 7 42 



laTA^e^ea;i° ( } , (2.21) 



2 [210 

where ctt is the Thomson scattering cross section and A^e is the electron number density. The expansion 
parameter 6^ is defined by 

9e^'^. (2.22) 
In deriving cqs. (2.12) (2.21), we have ignored 0{9^) contributions. Using eqs. (2.12) (2.21), one can 
show that the photon number is conserved order by order in terms of the expansion parameter 9e- 



We now apply the present result of the generalized Kompaneets equation to the Sunyaev-Zel'dovich 
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effect for clusters of galaxies. We assume the initial photon distribution of the CMB radiation to be 
Planckian with temperature Tq: 

"o(^) = ' (2.23) 

where 

Substituting eq. (2.23) and eqs. (2.12) - (2.21) into eq. (2.10), and assuming To/Te <C 1, one obtains the 
following expression for the fractional distortion of the photon spectrum: 

^t|? = ^-^^ + OeY, + 9lY, + elY, + elY,] , (2.25) 



Yo = -4+X, (2.26) 
Y, = + + S'(-^ + lx) , (2.27) 



2 5 10 V 5 5 

2 8 5 5 5 30 

. .'(-f.H,). 

15 2505 - 7098 --2 14253 -, 18594 -, 

= Y + ^^-^^ 

12059^, _ 128^ 
140 21 105 

-2/ 7098 14253 - 102267 -, 156767 -, 1216 -. 64-. 

54 ( 18594 205003 - 1920 1024 -3 

V 35 280 7 35 

. ^=(-li.?|.), 

135 30375- 62391-2 614727 -, 124389-. 

Y4 = H X X^ H X-^ X^ 

32 128 10 40 10 

355703 -5 16568 - g ^ 7516 - 7 22 - « 11 -9 

+ +To^^ -y^ +210^ 

^2 / 62391 614727 - 1368279 -2 4624139 -3 157396 -4 
-|- o I h -X — JC -h JC — JC 

V 20 20 20 80 7 

30064 -s 2717 -e 2761 

7 7 210 

^4 / 124389 6046951 - 248520 -2 481024 -3 15972 -4 
-|- o -|- JC — JC H JC — JC 

V 10 160 7 35 7 

18689 -5 



-9- 



where 



.6/ 70414 465992-, 11792 19778 ^3 



y = <JT j dtNe , (2.31) 

X = Xcoth^^y^ , (2.32) 

S ^ (2.33) 

sinh(^- 

Note that the analytic forms of Yq, Yi and Y2 in cqs. (2.26) (2.28) agree with the resuhs obtained by 
ChalUnor and Lasenby (1997). Finally we define the distortion of the spectral intensity as follows: 

An(X) 

Ai = ^ ; ; . 2.34 



3. ANALYSIS OF THE CONVERGENCE OF THE POWER SERIES 

We now carefully study the convergence of the analytic expressions of eqs. (2.25) and (2.34). In order 
to do the task, first of all, we integrate eq. (2.9) directly by numerical integration. We confirm that the 
total photon number is conserved with excellent accuracy (< 10^*^) in the numerical integration. We arc 
now ready to compare the present numerical results with those obtained by the analytic expressions of eqs. 
(2.25) and (2.34) for various X-Tg regions and investigate the accuracy of the analytic expressions. 



3.1. RAYLEIGH-JEANS REGION 

In the Rayleigh-Jeans limit where X — > 0, eq. (2.25) is further simplified: 

10 40 280 640 J ^ ' 

As is seen explicitly from eq. (3.1), the convergence of the power expansion is very fast in the X ^ limit. 
Furthermore, we show in Fig. 1 the T^-dependence of the spectral intensity distortion eq. (2.34) for X = 1. 
As is expected, the convergence is extremely fast for fcsTe < 50keV. Relativistic corrections higher than 
0{6l) terms are almost negligible in this region. So far the Sunyaev-Zel'dovich effects have been measured 
in the Rayleigh-Jeans region. Therefore one can reliably apply the analytic expressions of eqs. (2.25), (2.34) 
and (3.1) to the analysis of the observed data. 



An{X) 
no(X) 



-2yee 
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in passing we remark that the form of eq. (3.1) is meaningful only in an idealized situation. In order 
that the higher order terms have a physical meaning it is necessary that the electron distribution is rigorously 
given by the relativistic Maxwellian distribution eq. (2.6) with a precisely-determined temperature Tg. In 
the real observation, the electron temperature Tg has a significant amount of observational error. This 
thereby restricts the precision of the formula (3.1). 

3.2. X » 4 REGION 

As one can see from eq. (2.26), the leading order contribution Yq vanishes at X = 4 X. Therefore, 
higher order corrections become more important in this region. In Fig. 2 we have plotted the Tg-dependence 
of the fractional spectral distortion at X — 4. It is seen that the dash-dotted line is closest to the exact 
result. It should be emphasized here that the dash-dotted line is the contribution including only first four 
terms in eq. (2.25). The dashed curve which includes all terms in eq. (2.25) shows a poorer agreement 
with the exact result. This means that the power series expansion in 9e in this region is not convergent 
but asymptotic for large T^. We conclude that the analytic expression which includes up to 0{6l) terms is 
reliable for fcsTg < 15keV in the X » 4 region. We recommend that the analytic expression which includes 
up to O^Ys terms (dash-dotted curve) be used for the analysis of the observational data for 15keV < fc^Te < 
30keV, X « 4. 

3.3. WIEN REGION 

Wo now study the Wicn region whore X > A. As is mentioned earlier, the Sunyacv-Zcl'dovich effects 
have been so far studied obscrvationally in the A <C 1 region. However, the effects will be observed in the 
Wien region in the future. For an illustrative purpose, we show the Te-dependence of the spectral intensity 
distortion of eq. (2.34) at A = 8 in Fig. 3. The convergence is very slow. All curves are diverging quickly 
from the solid curve (exact result) for fc^Te > 30keV. We conclude that the analytic expression including 
up to 0{6^) terms is reliable for fcsTg < 15keV. In Figs. 4-5 we show A//y for UbT^ = 15keV and fc^Te = 
20keV, respectively. We confirm the good accuracy of the analytic expression for ksTg = 15keV. 



3.4. CROSSOVER FREQUENCY 



Finally we study the crossover frequency Xq, where the spectral intensity distortion vanishes. It 
is known that the accurate determination of the Xq values is extremely important for the study of the 
Sunyaev-Zel'dovich effects (Rephaeli 1995). In Fig. 6, we have plotted the Tg-dependence of Xq for 
ksTe < 50keV calculated by the analytic expressions and also by the numerical integration (solid curve). 
The numerical result is well approximated as a linear function of 6^ for ^B^e < 20keV. It starts to deviate 
from the linear form for ksTe > 20keV. We have fitted the numerical result as follows: 

Xo « 3.830 (1+ 1. 16746*6 - 0.85336*2) ^^ 2) 

The errors of this fitting function are less than 1 x 10~^ for < ksTe < 50keV. 

Rephaeli (1995) discusses the consequence of the rclativistic shift of the crossover frequency on the 
value of the peculiar velocity of the cluster measured with the use of the kinematic Sunyaev-Zel'dovich 
effect. For a cluster with ksTg = 13.8keV, one obtains an error of 20 km/s in the deduced value of the 
peculiar velocity for the SUZIE experiment corresponding to an accuracy of 1 x 10~^ in eq. (3.2), by 
making an interpolation of Rephaeli's estimation. Therefore, one concludes that the current level of the 
observational accuracy has not reached the accuracy provided by the theoretical formula (3.2). 

Another impliaction of the accuracy 1 x lO^'' of cq. (3.2) will bo the following. Let us consider a cluster 
with the plasma temperature ksTe — 10 kcV making a proper motion with the velocity v = 1000 km/s. 
Then the effect of this proper motion on the Sunyaev-Zel'dovich effect will be of order (v/c)''^ / (ksTc/mcP) 
= 5 X 10~^. Thus the kinematic effect of the proper motion of the cluster will not exceed the accuracy of 
eq. (3.2). 

Finally we estimate an error 6 of the power series approximation in the X — T^ parameter space. We 
define 

^=|l-7r^l, (3.3) 

where AI is given by cq. (2.34) with eqs. (2.25) (2.33) and A/„„„j is the result of the direct numerical 
integration of eq. (2.1). In Fig. 7, we have plotted the results for 6 — 1%, 5%, and 10%. It should be 
remarked that the sharp peaks in Fig. 7 correspond to the points of X where the analytic curve and the 
numerical curve are crossing. This can be easily seen in Fig. 5 that the dashed curve (analytic curve) and 
the solid curve (numerical curve) are crossing at these points of X. It is obvious from Fig. 7 that the power 
series approximation is extremely accurate up to 50keV for the Rayleigh-Jeans region {X < 1.5). The 
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approximation has a larger error in the X w 4 region where the leading order contribution vanishes. Finally 
the power series approximation is good up to 15keV for the Wien region. 

4. CONCLUDING REMARKS 

We have calculated the relativistic corrections to the Sunyaev-Zel'dovich effect including terms up to 
the order O{0l). The resuhs of the obtained analytic expressions have been compared with the result of 
the direct numerical integration. The extent of the applicability of these analytic expressions has been 
carefully examined. It has been shown that the present analytic expressions have an excellent accuracy 
in the Rayleigh- Jeans region where X < 1. There the convergence is extremely fast and the results are 
reliable at least for ksTe < 50keV. On the other hand, the applicability of the analytic expressions becomes 
to be limited for both cases of the X « 4 region and the Wien region where X > 4. In these regions, the 
presently obtained analytic expression which takes into account up to 0{9l) terms is reliably applicable 
for ksTe < 15keV. This is the highest electron temperature of the presently known clusters of galaxies. 
Therefore, the present analytic expressions can be reliably applied to the known clusters of galaxies. 

We thank Dr. A. Lasenby, Dr. A. Challinor, and Professor M. Longair for sending us a preprint and for 
stimulating discussions on this subject. We also thank our anonymous referee for many valuable comments 
which helped us tremendously in revising the manuscript. 
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Figure Captions 

• Fig.l. Spectral intensity distortion Al/y as a function of ksTe for X = 1. The solid curve shows 
the result of the numerical integration. The dotted curve shows the contribution of the first two 
terms in eq. (2.25). The long-short dashed curve shows the contribution of the first three terms. The 
dash-dotted curve shows the contribution from the first four terms. Finally the dashed curve shows 
the contribution from all the terms in eq. (2.25). 

• Fig.2. Same as for Fig. 1 but for X = 4. 

• Fig. 3. Same as for Fig. 1 but for X = 8. 

• Fig. 4. Spectral intensity distortion A//?/ as a function of X for fcsTe = 15keV. See Fig. 1 for the 
explanation of the curves. 

• Fig.5. Same as for Fig. 4 but for fc^Te = 20keV. 

• Fig. 6. Crossover frequency Xq as a function of ksTp. The solid curve shows the result of the 
numerical integration. The dotted curve shows the contribution of the first two terms in eq. (2.25). 
The long-short dashed curve shows the contribution of the first three terms. The dash-dotted curve 
shows the contribution from the first four terms. Finally the dashed curve shows the contribution 
from all the terms in eq. (2.25). 

• Fig. 7. Estimation of the error 5. The solid curve shows the results for 5 = 1%. The dash-dotted curve 

shows the results for 5 = 5%. Finally the dashed curve shows the results for 6 = 10%. 



